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Abstract 



We present total rates and kinematic distributions for the associated production of a single bottom quark and a Higgs 
boson at the Tevatron and the LHC. We include next-to-leading order QCD corrections and compare the results 
obtained in the four and five flavor number schemes for parton distribution functions. 



1 



Higgs boson production with one bottom quark jet at hadron colliders 

S. Dawson'^, C.B. Jackson'', L. Reina'', and D. Wackeroth'^ 
a) Department of Physics, Brookhaven National Laboratory, Upton, NY 11973-5000, USA 

b) Physics Department, Florida State University, Tallahassee, FL 32306-4350, USA 
c) Department of Physics, State University of New York, Buffalo, NY 14260-1500, USA 

We present total rates and kinematic distributions for the associated production of a single bottom 
quark and a Higgs boson at the Tevatron and the LHC. We include next-to-leading order QCD 
corrections and compare the results obtained in the four and five flavor number schemes for parton 
distribution functions. 



One of the most pressing problems of particle physics is 
to uncover the origin of electroweak symmetry breaking. 
The Standard Model (SM) predicts the existence of one 
scalar particle, the Higgs boson, as a consequence of the 
generation of gauge boson masses via the Higgs-Kibble 
mechanism. Extensions of the SM often introduce more 
Higgs fields, whose properties may drastically differ from 
the SM Higgs boson. Finding experimental evidence for 
the existence of one or more Higgs bosons and measuring 
their couplings to gauge bosons, leptons, and quarks is a 
major goal of particle physics. 

In this letter, we focus on Higgs boson production in 
association with bottom (6) quarks. The coupling of the 
Higgs boson to a 5 quark is suppressed in the SM by the 
small factor nib/v, where v— (\/2Gf)^^^^ = 246 GeV, im- 
plying that associated production of a SM Higgs boson 
with b quarks is very small at both the Tevatron and the 
LHC. In a two Higgs doublet model or a supersymmet- 
ric model, however, this coupling is proportional to the 
ratio of neutral Higgs boson vacuum expectation values, 
tan /?, and can be significantly enhanced for large values 
of tan /3. 

The associated production of a Higgs boson with b 
quarks proceeds at tree level via qq,gg — > bbh, as illus- 
trated in Fig. ^ Fully inclusive or semi-inclusive Higgs 
production cross sections are then obtained by imposing 
no identification cuts on the final state b quarks or by 
requiring that at least one b quark is observed at high 
transverse momentum, px, respectively. This approach 
with no b quarks in the initial state is dubbed the fixed or 
four flavor number scheme (4FNS). Large logarithms of 
the form A = log(^^/TO^) (for fih ~ M/i) arise at all orders 
in as from the integration over the px of the final state 
b quarks that originate from the coUinear splitting of an 
initial state gluon. In order to stabilize the perturbative 
expansion of the corresponding cross section, these log- 
arithms can be resummed into a b quark perturbatively 
defined parton distribution function (PDF), which inte- 
grates over the low px region of the corresponding b quark 
up to scales of the order of the factorization scale. This 
approach is identified as the variable or five flavor num- 
ber scheme (5FNS) 0,1111. The 5FNS is based on the 
approximation that the outgoing b quarks are at small 
Pt- It can be used only if the corresponding b quark is 



treated inclusively and no pr cut is required to identify 
the b quark in the final state. In the 5FNS the lead- 
ing processes for fully inclusive and semi-inclusive Higgs 
production with b quarks are bb h and bg —^ bh, re- 
spectively (the case of bg bh is illustrated in Fig. |21 at 
tree level). In this approach the processes qq,gg — > bbh 
are subleadingcontributions of higher order in the a" A"* 
expansion 0,13 ■ 

Assessing the validity and compatibility of the 4FNS 
and 5FNS approaches in the context of Higgs production 
with b quarks has recently been the subject of much theo- 
retical interest. In this letter, we focus exclusively on the 
production of a Higgs boson with one identified b quark 
jet, and analyze in detail the results obtained within the 
4FNS and 5FNS. Requiring one final state 6 quark mea- 
sures unambiguously the Yukawa coupling of the b quark 
and significantly enhances the rate with respect to the 
case when both b quarks are identified. Higgs boson pro- 
duction with one 6 quark jet (and h — > bb) has been 
extensively studied by the CDF and DO collaborations 
and is going to play a major role in the experi- 
mental searches for Higgs bosons beyond the SM at the 
Tevatron and at the LHC. Thus, a more dedicated effort 
aimed at refining the theoretical predictions for both to- 
tal and differential cross sections is mandatory. 

A first study of the next-to-leading order (NLO) QCD 
total cross sections for qq, gg — > bbh 0,0 and for bg — > bh 
pH l processes has been presented in Ref. • this let- 
ter we concentrate on the comparison of total and differ- 
ential cross sections at NLO QCD in the 4FNS and 5FNS 
schemes. This is the first comparison of differential cross 
sections in the two PDF schemes and is important to 
assess the residual theoretical uncertainties in future ex- 
perimental analyses. In particular, we discuss the effects 
of including the closed top quark loop diagram of Fig. |31 
a contribution that had been previously neglected, in the 
NLO calculation of bg bh in the 5FNS. 

The NLO QCD corrections to pp,pp b{b)h produc- 
tion in the 4FNS consist of calculating the 0{as) virtual 
and real QCD corrections to the qq, gg bbh tree level 
processes |^|^, imposing identification cuts on the trans- 
verse momentum and pseudorapidity of cither the 6 or 6 
final state quark (antiquark) . Results from the two ex- 
isting calculations 0, have been compared and found 
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FIG. 1: Sample Feynman diagrams for gg — > bbh and qq 
bbh production at tree level. 
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FIG. 2: Feynman diagrams for gb bh production at tree 
level. 



in good agreement (see Ref. Except for the iden- 

tification cuts, the calculation is identical to that for tth 
production El El [11 [H E [13 with the global re- 
placement of the top quark by the bottom quark mass 
{rrit — *■ rrif,) and Yukawa coupling {gtth ^ 9bbh)i where 
nih is always non zero. 

The NLO QCD corrections to pp,pp — » bh + bh pro- 
duction in the 5FNS have been presented in Ref. 1 1(1 and 
are encoded in the Monte Carlo program MCFM |18l|. In 
Ref. E3| > the calculation of the cross sections for bg bh 
is performed in the rub — approximation (except for the 
b quark Yukawa coupling), and for this reason the only 
virtual diagram containing a top quark loop (see Fig. I^)) 
is neglected. Indeed, the contribution of this diagram to 
the virtual cross section is proportional to gtth9bbh^b / "^it 
and therefore vanishes when = 0. At the same time, 
in the SM {gbih — i^b/v) the contribution of this diagram 
is overall of order g^j^^ as all other diagrams retained in 
the TO{, = approximation. So, it can play a relevant nu- 
merical role in the comparison between the 5FNS and the 
4FNS, where diagrams with closed top quark loops are 
included. To investigate this issue we have computed this 
contribution to the bg — > bh process and implemented it 

b b 



FIG. 3: Feynman diagram for the closed top quark loop con- 
tribution to gb bh. 
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FIG. 4: Total LO and NLO cross sections for pp,pp b{b)h 
production in the 4FNS as a function of /j, = fir ~ l-if for 
Mh = 120 GeV, at both the Tevatron and the LHC. 



into MCFM. All numerical results in the 5FNS presented 
here are obtained with this modified version of MCFM. 

Our LO numerical results are obtained using 
CTEQ6L1 PDFs 19, 20] and the 1-loop evolution of a^, 
while for NLO resuhs we use CTEQ6M PDFs and the 
2-loop evolution of as, with as{Mz) = 0.118. We use 
the A/5 running b quark mass in the b quark Yukawa 
coupling, evaluated at 1- and 2-loops respectively for LO 
and NLO results (with pole mass mi, = 4.62 GeV). Our 
renormalization scheme decouples the top quark from the 
running of mb{p) and agin) and is explained in detail 
in ReL We work in the SM but the results can be 
straightforwardly generalized to the case of the scalar 
Higgs bosons of a supersymmetric extension of the SM, 
for instance, by replacing the SM top and bottom quark 
Yukawa couplings accordingly, as discussed in Ref. . 

In order to simulate the experimental cuts, we require 
one of the final state b quarks to have px > 20 GeV and 
pseudo-rapidity 1 77 1 < 2.0 for the Tevatron and 1 77 1 < 2.5 for 
the LHC. In the NLO real gluon emission, the final state 
gluon and b quarks are considered as separate particles 
only if Ai? > 0.4 ( Ai? = ^/{KqfTiKW) ■ 

In Fig. ^ we show, for M;i = 120 GeV, the dependence 
of the LO and NLO total cross sections, calculated in 
the 4FNS, on the arbitrary renormalization/factorization 
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FIG. 5: Total NLO cross section for pp,pp — > b{b)h pro- 
duction at the Tevatron and the LHC as a function of Mh- 
We have assumed /ir =M/ = A'o/S for the central curves (see 
inlays) and varied fir and /i f independently to obtain the un- 
certainty bands, as explained in the text. The solid curves 
correspond to the 4FNS, the dashed curves to the 5FNS. 



scale fj, (with fir = H f = fJ-) ■ The NLO result has consid- 
erably less sensitivity to the scale choice, and the region 
around /iw/io/2 (/io = mf, + shows the least sensi- 

tivity to the variation of fi. For this reason we use /io/2 
as our reference scale in the following plots, whenever 
fir — Hf- Analogous results for the 5FNS total cross sec- 
tions have been presented in Ref. [loj . 

Fig. El shows the dependence of the NLO total cross 
sections on M^, in both the 4FNS and 5FNS. The bands 
illustrate the theoretical uncertainty due to the inde- 
pendent variation of fir and fif about the central value 
^r = /U/=A'o/2 (see inlays), between (0.2,0.25)^0 (4FNS, 
5FNS) and fiQ. It is extremely interesting to note that 
the 5FNS band is completely within the 4FNS band, and 
the corresponding central values are almost identical at 
the Tevatron and very close at the LHC. Including the 
closed top quark loop diagrams lowers the 5FNS cross 



section (by « 15% at the Tevatron and « 10% at the 
LHC, when Mh = 120 GeV and fir = = 0.5/io) and 
makes the theoretical prediction in the 4FNS and 5FNS 
fully compatible (see for comparison Fig. 6 in Ref. [lll|). 
Note that the bands only give an indication of the theo- 
retical uncertainty of each approach due to the residual 
scale dependence. Other sources of theoretical uncertain- 
ties, like PDF uncertainties, have not been considered. 
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FIG. 6: da/dpT at the Tevatron and the LHC for Mh = 
120 GeV and fir — fif = /io/2. We show the NLO results in 
the 4FNS (solid) and 5FNS (dashed), using two different bin 
sizes, 2 GeV and 12 GeV (see inlay). 



Finally, in Figs. I6I8I we compare the results for the 
Pt and rf distributions of the Higgs boson in both the 
4FNS and 5FNS, at the Tevatron and the LHC. We see, 
in general, a good agreement between the two schemes, 
except in regions of kinematical boundaries. This is par- 
ticularly dramatic in the distributions where, around 
~ 20 GeV, a kinematical threshold causes the 5FNS 
NLO calculation to be highly unstable. This instability 
can be reabsorbed by using a larger bin size (see inlays) , 
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and could therefore be interpreted as a sort of theoret- 
ical resolution for the 5FNS. The instabilities could be 
removed by a systematic resummation of threshold cor- 
rections '23,113, but this is not implemented in MCFM. 
Fig. |S1 illustrates the impact of NLO QCD corrections on 
p'^ and rjfi distributions in terms of a differential K-factor 
[da N LO / d(T Lo) ■ It is interesting to note that the 4FNS 
and 5FNS agree at large but they differ substantially 
at low py. As can be seen in Fig.|Hl there are regions of 
and rjh where the NLO QCD corrections can consid- 
erably affect the shape of the distributions. 

In this letter we have shown how the theoretical pre- 
dictions for both total and differential cross sections for 
pp,pp b{b)h production within a 4FNS and 5FNS are 
fully compatible within the existing theoretical errors due 
to the residual normalization and factorization scale de- 
pendence at NLO in QCD. This is crucial to experimental 
searches based on bbh production when only one b quark 
jet is identified. 




FIG. 7: da/dr/h at the Tevatron and the LHC for Mh = 
120 GeV and fir = l-if ~ Mo/S. We show the NLO results in 
the 4FNS (solid) and 5FNS (dashed). 
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